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In  this  study,  perturbation  method  is  applied  to  obtain  the  solution  of 
dynamical  system  which  does  not  involve  periodic  solution.  Since 
perturbation  method  presents  results  in  series  approximation  form, 
accuracies  with  retention  of  different  number  of  terms  are  studied. 
Likewise,  different  nonlinearities  are  also  considered,  with  the  level 
of  agreements  of  their  solutions  are  also  presented  in  tabular  form  to 
facilitate  ease  of  comparison.  Numerical  results  obtained  in  all  cases 
revealed  that  retention  of  one  or  two  terms  in  the  solution  suffices  for 
a  reasonable  level  of  accuracy  for  system  with  no  periodic  solution. 
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1.  Introduction 

Many  dynamical  systems,  when  modeled  result  into  differential  equations,  mostly  nonlinear 
ordinary  differential  equations  [1,2].  Dynamical  system  such  as  simple  pendulum,  conservative  and 
non-conservative  oscillators,  and  so  on,  belong  to  the  dynamical  systems  with  periodic  solutions. 
Many  researchers  have  worked  on  the  solution  of  this  category  of  problems  and  proffer  reliable 
solutions  to  both  conservative  and  non-conservative  oscillators.  Prominent  among  the  methods  of 
solution  are;  modified  differential  transform  method  for  solution  of  excited  nonlinear  oscillators 
under  damping  effects  by  [3],  nonlinear  oscillator  with  discontinuity  was  discussed  by  [4],  where 
He’s  energy  balance  method  was  adopted,  [5]  used  analytical  approximation  technique  to  solve 
strongly  nonlinear  oscillator  problems,  just  to  mention  a  few. 

While  on  the  other  hand,  there  is  another  category  of  dynamical  systems  that  do  not  involve  periodic 
solutions.  Such  problems  work  mainly  on  the  principles  of  moment,  where  loads  are  placed  on  static 
beams  fixed  at  both  ends,  is  a  good  example  of  dynamical  system  that  does  not  involve  periodic 
solution.  A  beam  is  a  structural  item  with  a  length  that  is  considerably  longer  than  its  width  and 
thickness  [6].  Solution  to  such  problems  were  considered  in  [1,6,7, 8].  While  perturbation  and 
modified  perturbation  methods  were  employed  in  [6]  to  solve  dynamical  problems  with  no  periodic 
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solution  and  nonlinear  oscillators  respectively,  [3]  on  the  other  hand  used  modified  perturbation 
theory  for  anharmonic  oscillator. 

In  the  present  work,  the  approach  proposed  in  [7]  is  adopted  while  using  perturbation  theory  to 
handle  dynamical  problems  with  no  periodic  solution,  but  with  varying  degrees  of  nonlinearities.  It 
has  been  shown  here  too,  that  after  injecting  the  results  of  asymptotic  expansion  in  the  given 
problem  and  resulting  system  of  nonlinear  initial  value  problems  solved,  retention  of  few  terms  in 
the  final  series  solution  suffices  for  an  appreciable  level  of  accuracy. 


2.  Methodology 

2.1.  Perturbation  Method 

Perturbation  Method  (PM)  involves  the  introduction  of  a  perturbation  term  into  the  nonlinear  term 
that  occurs  in  the  given  differential  equation,  if  it  does  not  come  with  the  problem.  The  perturbation 
term,  which  depends  on  a  small  parameter,  s,  that  is  added  to  the  system  to  bring  about  corrections 
[6].  These  corrections  which  are  mostly  smaller  compared  to  the  size  of  the  other  quantities,  are 
now  calculated  using  asymptotic  series.  This  method  gives  best  results  when  applied  to  non- 
dynamical  system  problems  with  a  non-periodic  solution  (static  systems).  To  illustrate  the  method, 
we  shall  use  the  nonlinear  initial  value  problem 

y"(x)  -f  y(x)  =  y(0)  =  A,  y'(0)  =  0  (1) 

To  solve  (1)  by  PM,  a  small  nonnegative  perturbation  term  e  is  introduced  to  the  nonlinear  term 
y'(x)^.  Thus,  the  IVP  in  (1)  becomes 

y''(x)  -f  y(x)  =  £y'(x)^  y(0)  =  A,  y'(0)  =  0  (2) 

The  asymptotic  expansion  of  y(x,  s)  is  given  as 

n 

y(x,£)  =  ^  £''y/c(x)  +  0(£”+^)  (3) 

k=0 


The  equivalent  form  of  (3)  is 

y(x,  £)  =  yo{x)  +  £yi(x)  +  £^y2(x)  -I-  — h  e^y-^ix)  +  (4) 

Either  (3)  or  (4)  is  substituted  into  (2),  to  get 

y'oix)  +  eyiix)  +  s'^y 2  {x)  +  0{£^) 

=  ^(yoW  +  £y{  W  +  £^yK^)  +  0(£^))  -  yo(x)  -  £yi(x)  -  £^y2(x) 

+  (5) 

y'oix)  +  £y['(x)  +  £^y2  W  +  •■• 

=  ^(yoW^  +  2£yo(Y)yj(x)  +  2£2yi'(x)2  +  2£Vo'Wy2'W  +  ■■•)  -yoW 
-£yi(x) -£2y2(x)  ...  (6) 
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yo  (y)  +  £y['(x)  +  £^y2  W  +  •■• 

=  -yo W  -  £(yi W  -  yo'W^)  -  £^(y2 W  -  2yo(x)y;(x)) 

-  £^(y3  W  -  2yo(x)yKx)  -  yi'(x)^)  +  (7) 

Comparing  the  coefficients  of  various  power  £,  the  IVPs  are  obtained: 

y^'(x)  +yo(x)  =  0,  y(0)  =  A,  y'(0)  =  0  (8) 

£^:  y['(x)  +yi(x)  =  yo'(x)^  yi(0)  =  0,  y/CO)  =  0  (9) 

yz'Cx) +y2(x)  =  2y^(x)y;(x),  y2(0)  =  0,  y'(0)  =  0  (10) 

y^'(x)  +  ygCx)  =  2y^(x)y^(x)  +  yi'(x)^  ysCO)  =  0,  y^(0)  =  0  (11) 

What  follows  is  the  solution  of  the  system  of  Equations  (8),  (9),  (10)  and  (11)  using  any  of  the 
known  methods.  We  therefore  obtain  yo(x),yi(x),y2(x),  ...  and  substitute  them  in  (4)  to  get  y(x) 
as 

y(x)  =  yo(x)  +  £yi(x)  +  £^y2(x)  +  (12) 

Sometimes,  solving  for  yo(x),yi(x)  and  y2(x)  suffices,  and  some  other  times,  solving  for  the  first 


two  may  suffice.  As  stated  earlier,  if  the  problem  being  solved  is  a  dynamical  system  that  does  not 
involve  periodic  solutions  then  solving  for  yo(x)  alone  or  at  most  solving  for  yo(x)  and  yi(x) 
suffices. 

3.  Numerical  Experiments 

To  illustrate  the  implementation  of  Perturbation  Method  (PM)  described  in  the  earlier  sections  of 
this  paper,  we  consider  a  static  system  that  does  not  have  periodic  solutions,  since  this  where  PM 
gives  optima  results,  [6]. 

Probleml  [6] 

Consider  the  structural  beam  with  two  hinged  ends  under  a  concentrated  load  at  the  middle  of  the 
span.  The  governing  equation  for  this  system  is  given  by  the  following  IVP: 

Fx  3 

y"(x)  “  +  y'W^)^  =  O'  (13) 

together  with  the  initial  conditions 

y(0)  =  0,  y'(0)  =  (14) 

where  y(x)  is  the  beam  deflection  (vertical  displacement),  y'(x)  is  the  corresponding  slope  and  x 
is  the  time.  The  constants  E  and  I  are  the  modulus  of  elasticity  and  the  moment  of  inertia  of  the 
cross  section  about  its  axis,  respectively.  E  is  the  length  of  the  beam  and  E  is  the  load  which  is 
concentrated  at  the  middle  of  the  span. 

Solution 

In  this  case  the  problem  did  not  come  with  the  required  small  parameter,  so  it  shall  be  introduced 
by  adding  it  to  the  nonlinear  term  as  follows: 

Fx  s 

y"(x)  “  ^  (1  +  £y'(^)^)2  =  0,  (15) 

FI? 

With  the  conditions:  y(0)  =  0  and  y'(0)  =  —  using  the  expansion: 

2.E1 

n 

y(x,  f)  =  ^  £^y/c(x)  +  0(f”'''^)  (16) 

k=0 
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In  (3),  when  (3)  is  written  in  the  form 


(y"W)^=  ^(1+ 

We  have 


(Y£'^yk(x)  +  0(s^+^)\  =  1+  W + 

Vk=0  /  Vk=o  / 

[yo  W  +  ^yl'W  +/y2ix)  +  o(8^)r 

=  [1  +  £(yoW  +  £yi  W  +  £^yK^)  +  ■■•  YY 


(yo  W)^  +  2£yo  Wyi  W  +  f^yi  +  2yo  Wy2  W)  +  •■• 

=  +  3£(y^(x))^  +  £2  (3(y^(x))'^  +  6y^(x)y;(x)  j 


((yoW)^  +  12(y^(x))V'iW  +  3(y{(x))^  +  6y^(x)y^(x)j  +  •■•  (20) 


Equating  the  eoeffieient  of  £°,  and  in  (8)  gives: 


2  FLr 

£“:(y;w)  =  JPT?.  y(0)  =  0,  y;(0)=^ 

3F^x^  7 

£y  2y^'(^)yi'W  =  4^2/2  (yoW)  <  yi(o)  =  o<  y'o  =  o 


£2;  (y"(Y))^  +  2y^'(x)y^'(^)  =  ^^{^iyo(^)T  +  2yo(Y)yi(x) j ,y2(0)  =  y'^CO)  =  0  (23) 


Solving  (21)  and  (22)  shall  suffiee,  as  the  solution  to  the  (23)  will  be  eumbersome  and  may 
therefore  not  be  neeessary. 

From  (21), 

Fx  FL^ 

yo  W  =  ^  .  yo(0)  -  0,  y'„(0)  =  ^  (24) 


y  \j  2iEI  ^  '  y  o 

Fx^ 

y;w  =  i^+‘'i’ 

Using  the  assoeiated  eondition,  we  have 


Substituting  for  in  (25),  we  get 
Fx^  FL^ 
yo(^^  “  4£7  ^  16EI 

Fx^  FL^x 
~  12EI  ^  16EI  ^ 


Using  the  eondition  yo  =  0,  we  have 

C2  =  0. 
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Thus,  (27)  becomes 

Fx^  FL^x 

“  125/  ^  165/ 


(28) 


(28)  is  taken  as  the  solution  of  (21),  not  losing  focus  on  the  fact  that  (21)  will  have  two  answers 
that  only  differs  in  sign.  Since  (28)  satisfies  the  subsidiary  conditions  in  (21),  it  is  accepted  as  the 
solution. 

We  now  proceed  with  the  solution  of  (10)  as  follows: 

3F^x^  2 

yi  W  =  ^^(y'o(22))  ,  yi(0)  =  y\(0)  =  0 


yi(x) 


3F^x^ 

4^2/2 


3F^x^  (Fx^  F\}  \ 
■  452/2  +  165/ ) 

FI  (Fx^  ^  FL^  \ 
V4£7  ^  165/y 


yi'(x) 


3Fx  fFx^  ^  FL^  \ 
4£7  \4K[  ^  165/ ) 


3F^x'^  3F^L^x^ 

~  6452/2  +  12852/2  +  ^3 

Using  the  associated  conditions,  (22)  gives 

C3  =  0. 

Thus, 

^  _  35^x4  ^  35^12x2 
~  6452/2  +  12852/2 

And 

352^5  F^L^x^ 

~  32052/2  +  12852/2  +  ^4 

Using  yi(0)  =  0,  (30)  gives 

C4  =  0. 

Thus, 

3F'^x^  3F^L^x^ 

~  32052/2  +  12852/2 

Hence,  the  solution  of  (15)  is  given  by 

yW  =  yoW  +  £yiW  +  •■• 

As 

_FL^x  Fx^  /3F^L^x^  3F^x^  \ 

~  165/  ^  125/  ^  ^  1^12852/2  ^  32052/2  J  + 


(29) 


(30) 


(31) 


(32) 


(33) 


Since  £  is  a  very  small  nonnegative  constant,  taking  £  =  1  in  this  case  yields 
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y(x)  = 


FL^x  ^  Fx^  ^  3F^L^x^  ^  3F^x^ 

16E1  ^  125/  ^  128^2/2  +  320^2/2  ■ 


Problem  2  [6] 

Solve  the  BVP  in  Problem  1  by  setting 


(1  +  y'(x)^)2  X  1  +  2 


And  compare  your  results  with  these  obtained  in  Problem  1. 


(34) 


Solution 

The  new  problem  is 

(l  +^y'(x)2)  =  0, 


y”(x) 


Fx 


2E1 

with  the  conditions 
y(0)  =  0  and  y'(0)  = 


FL^ 

16E1 


(35) 


(36) 


To  solve  this  problem  by  permutation  method,  the  small  positive  parameter  s  is  introduced  to  the 
nonlinear  term  in  (35).  Thus,  we  have 


Fx  /  3 

y"(x)-  — (1  + 


2EI 


-fy'(x)^)  =  0 

Using 

n 

y(x,  f)  =  ^  £^yfc(x)  -I- 


k=0 

in  (35),  we  have 


n  /  n  \  ^ 

^  fVk  (x)  +  0(£”+^)  1 +  0(£”+^)  J 

k=0  V/c=0  / 


(37) 

(38) 


(39) 


To  (x)  +  £yi  (x)  +  £^y2'(x)  +  0(£) 


Fx 


1  +  2  ^(yo(^)  +  ^yi(^)  +  ^^yK^)  +  o{sfY 


yo 


(x)  +  £y{'(x)  +  £^y2'(x)  +  ^  +  3£(yo^(x))  ^ 


+ 


Taking  the  coefficients  of  various  powers  of  £,  we  have 

Fx  EL? 

y"oM  =  j^.  yo(0)  =  0,  y;(0)  =  j^ 
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3Fx 

£1;  y'i(x)  =  yi(o)  =  0.  yi(o)  =  o 

3Fx 

s^:  y'^(x)  =  :^yo(Y)y{(x),  y2(0)  =  0,  y^(0)  =  0 

2E1 

Solving  (42),  (43)  and  (44)  successively,  we  have 

Fx 

YoM  = 


2EI 


y'oW 


Ex^ 

4^ 


+  C-i 


Using  the  associated  eonditions,  we  get 

EL^ 

Cl  = 

Hence  (46)  becomes 


16EI 


Ex^  EL^ 

~  4£7  ^  16EI 

Ex^  EL^x 
~  12EI  ^  16EI 


+  C2 


Using  yo(0)  =  0,  we  have 


0, 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


hence  (49)  becomes 

Ex^  EL?x 
“  125/  ^  165/ ■ 

For  the  solution  of  (43),  we  proceed  as  follows: 


y'i(x)  = 
y'i(x)  = 


3Ex  fEx^  ^  EL^\ 


4EI  \4EI 

3E^x^ 
16F2/2  ^ 


16F/, 


3E^x^ 

y;cx)  =  g^+ 

Using  the  attached  eonditions,  we  have 

C3  =  0. 


(50) 


3E^x^L^ 

^  SF^L^x 

128F3/3 

^  64F3/3 

(51) 

3E^x^L^ 

^  3E^L^x^ 

S12E^P 

^  128F3/3  ^ 

(52) 

Thus  (52)  reduces  to 
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_  3F^x^  SF^L^x^ 

~  64^2/2  +  512^3/3  +  128^3/3 

SF^x^  ^  SF^x^L^  ^  F^L^x^ 

~  320F2/2  ^  3060F3/3  ^  128E^P  ^ 

Using  the  associated  condition  here  again  gives 

C4  =  0. 


Hence, 


3F2x^  ^  SF^x^L^  ^  F^L^x^ 

~  320F2/2  ^  3060F3/3  ^  128E^P 


(53) 

(54) 


(43) 


Solving  (44)  becomes  cumbersome,  hence  the  first  two  results  shall  suffice  here  again.  Therefore, 


yW  =yoW  +  £yiW  +  •■• 


y(x)  = 


Fx^  FL^x 

- 1 - h  £ 

12EI  16E1 


(  3F^x^  3F^x^L^  F^L^x^  \ 

l^320F2/2  ^  3060F3/3  ^  128E^P ) 


+  •■• 


Using  f  =  1,  we  have 

,  ^  _  Fx^  ,  ,  ^F^x^  ^  3F^x^L^  ^  F^L^x^ 

~  12EI  ^  16EI  ^  320F2/2  ^  3060E^P  ^  128E^P' 


(44) 

(45) 


(38) 


4.  Results  and  Discussion 


Table  1:  Results  for  y(x)  =  yo(x)  for  Problems  1  and  2 


X 

L 

El 

F 

Problem  1 

Problem  2 

1 

1 

500 

100 

0.029166710 

0.029166710 

2 

1 

500 

100 

0.15833310 

0.158333100 

3 

1 

500 

100 

0.48751000 

0.48751000 

4 

1 

500 

100 

1.11667100 

1.11667100 

5 

1 

500 

100 

2.145583100 

2.14583100 

Table  2:  Results  for  y(x) 

=  To  (^)  +  Ti  (^)  for  Problems  1  and  2 

X 

L 

El 

F 

Problem  1 

Problem  2 

1 

1 

500 

100 

0.029854210 

0.02961210 

2 

1 

500 

100 

0.172833100 

0.17108410 

3 

1 

500 

100 

0.587063100 

0.58221810 

4 

1 

500 

100 

1.52067100 

1.51271000 

5 

1 

500 

100 

3.35677100 

3.35003100 

Table  3:  Results  for  y(x) 

—  ToCy)  for  Problems  1  and  2 

X 

L  El 

F 

Problem  1 

Problem  2 

1.0 

3  500 

100 

0.12916710 

0.12916710 

2.0 

3  500 

100 

0.35833310 

0.35833310 

115 


B.M  Yisa  /  NIPES  Journal  of  Science  and  Technology  Research 


2(2)  2020  pp.  108-117 

3.0 

3 

500 

100 

0.7875100 

0.78751000 

4.0 

3 

500 

100 

1.51667100 

1.51667100 

5.0 

3 

500 

100 

2.64583100 

2.64583100 

Table  4: 

Results  for  y(x) 

=  yo(^)  +  yi(^)  for  Problems 

1  and  2 

X 

L 

El 

F 

Problem  1 

Problem  2 

1 

3 

500 

100 

0.13235410 

0.1313100 

2 

3 

500 

100 

0.39283310 

0.38609210 

3 

3 

500 

100 

0.95456310 

0.94134100 

4 

3 

500 

100 

2.08067100 

2.08095100 

5s 

3 

500 

100 

4.16927100 

4.24923100 

Table  5: 

Results  for  y(x) 

=  yo(x)  for  Problems  1  and  2 

X 

L 

El 

F 

Problem  1 

Problem  2 

1 

5 

500 

100 

0.32916710 

0.32916710 

2 

5 

500 

100 

0.75833310 

0.75833310 

3 

5 

500 

100 

1.38751000 

1.38751000 

4 

5 

500 

100 

2.31667100 

2.316671000 

5 

5 

500 

100 

3.6458310 

3.64583100 

Table  6: 

Results  for  y(x) 

=  yo(^)  +  yi(^)  for  Problems 

1  and  2 

X 

L 

El 

F 

Problem  1 

Problem  2 

1 

5 

500 

100 

0.33735410 

0.33755100 

2 

5 

500 

100 

0.83283310 

0.83910810 

3 

5 

500 

100 

1.68956100 

1.7372110 

4 

5 

500 

100 

3.20067100 

3.40145100 

5 

5 

500 

100 

5.79427100 

5.40702100 

4.1.  Discussion 

The  two  problems  considered  are  different  only  in  the  nonlinearity  parts  as  indicated  in  Problem  2. 
The  problems  were  solved  by  perturbation  method  with  the  usual  introduction  of  the  perturbation 
parameter,  s,  to  nonlinear  terms.  The  solutions  are  presented  by  retaining  the  first  term  if  the  series 
only  for  the  two  problems,  and  also  by  retaining  the  first  two  terms  of  the  series.  Solving  for  the 
third  term  in  the  two  cases  presented  very  cumbersome  situations  and  were  therefore  abandoned. 
Tables  1  and  2,  3  and  4,  and  5  and  6  presents  the  numerical  results  obtained  using  the  same  set  of 
values  assigned  to  the  terms  and  L  taken  three  values  respectively.  All  the  values  assigned  to  the 
constants  are  adopted  from  [6].  It  is  easily  seen  in  the  tables  that  the  solutions  with  only  yo(Y) 
retained  produced  the  same  results  for  the  two  problems,  as  we  have  them  in  Tables  1,  3,  and  5.  On 
the  other  hand,  with  yo(Y)  and  yi(x)  retained  the  results  are  still  close,  but  not  completely  the  same. 
As  stated  earlier,  retaining  one  or  two  terms  in  the  series  suffices,  since  there  are  no  significant 
differences  in  the  results  with  one  term  and  two  terms  retained. 

5.  Conclusion 

Perturbation  Method  (PM)  has  been  shown  to  have  performed  very  well  for  the  problems  considered 
in  this  work.  The  results  obtained  for  the  two  nonlinearities  are  very  close  and  satisfactory  results  are 
realized  by  retaining  one  or  two  terms  in  the  solution  series.  It  is  therefore  sufficient  to  affirm  that 
PM  is  very  suitable  for  solutions  of  dynamical  systems  that  do  not  possess  periodic  solution. 
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